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Introduction 


The pinhole-occulter system is a Space Shuttle based experiment for 
the production of hard X*ray images taken primarily fnmi the Sun. The 
systm is basically a pinhole cunera utilizing a deployable S0*m flexible 
boom for separating the pinhole ft^ the recording devices located In the 
Shuttle as seen In Figure 1. At the distal end of the boom from the 
Shuttle is a 50 kg mask containing pinholes and ccronograph shields. At 
the proximal end the detectors are located and munted, along with the 
tK)om, to a gimbal pointing system (either IPS or AGS) aligned with the 

A 

target along vector d. 

A 

The mask must be pointed at the X-ray source, along vector m with 
a iiigh degree of pointing stab'^llty to align the axes of the detectors 
with the pinholes and shields. Failure to do so will result In a blur- 
ring of the Images on the detectors and a loss of resolution. Being a 
Shuttle based experiment, the system will be subjected to the distur- 
bances of the Shuttle. The worst of these Is thruster firing for orbit 
correction; the Shuttle uses a bang-bang thruster control system to 
maintain orbit to within + 0.1®. Other disturbances Include man motion, 
motion Induced by other systems, and gravity gradient torques. 

The control system of the pointing mount can sense both position and 
velocity of the mask tip and uses these to accurately estimate the flex- 
ible body modes of the system. An optimal control /suppression scheme Is 
then used to control these nwdes. Disturbances are detected by 
sensors and are used as commands to drive the system. The AGS, with 
perfect sensors was utilized In this three axis control system study. 

The analysis of sensor/drive errors as well as a free body analysis will 
be the subjects of a later report. 
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i. 


FIGURE 1 
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Background 


The basis of the dynamic model was to usi* standard modal analysis. 
On any flexible structure, modes (eigenvalues) are set up and damp out. 
Each mode (a pair of complex conjugate eigenvalues) has two complex con- 
jugate eigenvectors, associated with It. By using a finite e1«nent model 
of the beam, there Is a mass associated with each node (connection points 
of the e1«nents). The eigenvectors along with these nodal masses can be 
used to calculate the deflection of the beam at any node due to any mode. 
Motions of the base excite th.' ''.odes and Induce deflections of the beam 
tip. 


Given these flexible body modes, the state equations for the system 
can easily be written as: [1,2] 

x • Ax + Bu (1 ) 

y « Cx + Du (2) 


The inputs to the system are accelerations 1n each of the three axes 
while displacements feed forward to the output. Partitioning the states 
of the system Into controlled states, x^ and suppressed states, x^, we 
have : 



(3) 




+ Du 


( 4 ) 
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By noticing that the double Integration of the acceleration Inputs 
are part of the outputs, we can define six (6) new states, Xj^, which are 
the rigid body trades of the boom. [3] This fact Introduced into Eqns. 3 
and 4 yields the equations of state for the system: 



Where the control vector u and output y are 





(5) 


( 6 ) 


three dimensional vectors. 


This parti on ing of the state vectors reflects the objectives of the 
control system design. The .controlled states, x^, must be controlled to 
achieve satisfactory system performance while the suppressed states, x^, 
are known but not critical to the control design [4,5]. 

What we would like to do is effect changes in x^ and no changes in 
either x^ or Xj^. The effect of the controller on the rigid modes or 
suppressed modes is called control spillover. We wish on one hand to 
optimize the System for the controlled states and limit the control 
spillover: the controlled states are controlled while the suppressed 
states are not excited by the controller and the rigid states are 
unaffected by It. 
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Let us procted as follows: suppose we have an optiiwm controller 


11 ° • - KXg (7) 

determined by using an Index of performance 

PI • ’ (xj Qj Xj ♦ Rj u) dt . (8) 

0 

We would then have In the control equations: 



Clearly, we would like to have for no control spillover: 


BgK • 0 

(10) 

Bj^K • 0 

(11) 


■.owever, equation *n implies K » 0 because of the form of B,^C3]. Equation 
10 implies that and K are orthogonal [5,6,7], VJe can modify equations 
10 and 11 so that 

K-0 , (12) 

which merely implies the orthogonality of B^ and K. 

By modifying the performance index to include spillover terms (B^ii) 
we can simultaneously minimize the ?I and the product B^K. Including 
these spillover terms in the PI, we have: 



( 13 ) 


PI ■ / [xj Qg Xg ii + (3^u)‘ Rj {3j u)] dt 

0 

■ * LxJ Qj X + u''' (Rj + S,' S, S, ) u] dt . (H) 

0 

Now by heavily penalizing the PI to suppressed state spillover (R^ »>Rg), 
we can force the maximization of the orthogonality of and K. The feed- 
back coefficient matrix can now be solved from the familiar matrix riccati 
equation of Appendix 2. 

The control spillover into the rigid body states cannot be avoided, 
but represents only steady state following errors. These errors can be 
lessened by the proper selection of controller gains for connand inputs. 

By the design of an optimal tracker, ort.hcganal to the suppressed states, 
the steady state following errors can be minimized. 

Once the optimal controller is determined, the boom tip responses 
and control inputs {-Kx ) can be calculated using the system equations 5, 

6 and 9 . Taking the LaPlace transform of equation 9 substituted into 


equation 3 we have: 


sXj(s) • (A^ - BjK) Xjts) + U(s) 

(15) 

5 X 5 ( 5 ) - (Aj) X5(5) - BjK Xj(s) * 5^ 85 U(5) 

(15) 

sX|,(s) - - B„K X|.( 5 ) + A„ X„( 5 ) + 5 ^ 8 ^ U(s) . 

(17) 

Rearranging the terms of equations 15-17 and defining 


Oq(s) • (SI - Aj. B^K)"'’ 


?s(s) • (SI - Aj)’^ 


5r(s) • (SI - A^)"'' , 

(13) 


we have: 
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(19) 


Xj($) ■ >,($) fj U(s) 

Xj(i) • «,(s) B, CI3 - X >o(») 5j] U($) (20) 

X,($) ■ »„(*) B„ CI 3 - X },($) Bj] U(»). (21) 

The term. K «q(s) In equations 20 and 21 represent, the control spill- 
over In the suppressed states and rigid body states respectively. This 
term represents the activation of the suppressed states in equation 20 
and a reduction of the magnitude of t.he response of the controlled states. 
In equation 21, this term represents s further reduction in the magnitude 
of the controlled states and the Introduction of steady state following 
errors. 

Using the output equation 6 for t.he syste.m, we can find the boom 
tip responses as: 


Y(s) 


C^.Xc(s) + CjXjls) ^ C;^X.(s) 


But 


K^qB^] )U(s). 
( 22 ) 


s^C|^S|^Bj^ » Ij since the positional Input Is fed forward to the 
output In eqn. 2. 

Therefore, 

ns) • Cs^C^S^B, + s^CjSjBjCj - K!„B,) . 13 - KS,B^]U(s) . (23) 

which are the controlled modes, supp'-essed -odes with spillover and 
rigid body modes with spillover. 
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The Tracking Problaw 

The tracking capability of the system represents the ability of the 
system to follow steady state conmand signals. In the case of the POF 
boom It Is desired to track second (2nd) order polynomials. T3] The 
systen out|HJt (tip angle of the boom) Is given by equation 23 but the 
only steady state terms are given by the rigid body iiwdes with spillover. 
Namely. 

Yjs^s) ■ (I 3 ‘ ^24) 


where U(s) Is the coninand signal. 

It Is desired to track a system of constant equation R. such that 


Y(s) 


R 


(25) 


with a command signal 
U(s) 


SjS + 02S + 


(26) 


which yields the tracking error 

2 

•V a<jS + a^S + 31 B 

E(s) • Y„ . Y(s) • (I 3 - ^ -) - ^ . (27) 

The first term of equation 27 reduces to a ratio of polynomials In s, 
whose denominator Is cwnposed of the optimum controlled modes. D(s). The 
total tracking error Is then the ratio of two polynomials In s; 


E(s) 


HiiL 

s^D(s) 



K, Ni(s) 
JL 4. — 

s 0 (s) 


(28) 


where the right half equation Is the partial fraction expansion of E(s). 

Notice that the steady state tracking errors are given by the first three 

N,(s) 

terms while the term, -p^ . represents the excited controlled modes 
which rapidly damp out. 
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The approach taken In this work was to varying the apa2 ”3 
matrices to obtain zero steady state error coefficient; namely Kq * ■ 

K2 ” 0. A complete derivation of this approach Is sN)wn In appendix one. 


Results 

The matrix riccatl equation was solved using the PI of equation 14 
and no control spillover occurred In the Z axis. This yielded an equation 
24. 

r 

1 - 


Y„(s) 


1 


1.765+2.5S,. . 
.99 

* 1 - 


0 

0 


4. 3+1. 2s 

. 0 

* ^ 

s^+2.05s+56.3 

1.765+2. 66s 

4. 3+1. 3s 

. 0 

s^+2.69s+2.99 

s^+2.05s+56.3 

0 

1 

1 

. 1 

found the tracking error could 

be 


(29) 


minimized by: 


and* 


U*(s) ■ 

X' ^ 13 xo 


,/.i - .828s^+3.408s+3 

,($) -3 


Uj(s) . -f . 


(30) 

(31) 

(32) 


Plots of e'(t), and e*(t) are shown In Figure Two. Following errors start 
^ y 

at zero as the shuttle starts drifting through Its deadband and gradually 
Increase to atout .1 arc sec after 100 sec. During a typical shuttle 

A 

deadband drift lasting 50 secs the maxlmun tracking error In x Is .04 arc 
sec and .03 arc sec In y. 
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inclusions 


Tht ginbtl driven R)F bo<»i his been shown to be capable of track- 
ing during siMttle gravity gradient drifting to within a few hundredths 

m 

of an arc sec. The regulator system was reoptimized to eliminate z 
axis tracking errors entirely. This reoptimization did not change the 
characteristic eigenvalues of the system. 
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APPENDIX I: DERIVATION OF TRACKING CONTROLLER 
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